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ABSTRACT 

We consider the possibility that the excess of cosmic rays near ~ 10 18 eV, reported by 
the AGASA and SUGAR groups from the direction of the Galactic Center, is caused 
by a young, very fast pulsar in the high density medium. The pulsar accelerates iron 
nuclei to energies ~ 10 20 eV, as postulated by the Galactic models for the origin of 
the highest energy cosmic rays. The iron nuclei, after about 1 yr since pulsar for¬ 
mation, leave the supernova envelope without energy losses and diffuse through the 
dense central region of the Galaxy. Some of them collide with the background matter 
creating neutrons (from desintegration of Fe), neutrinos and gamma-rays (in inelas¬ 
tic collisions). We suggest that neutrons produced at a specific time after the pulsar 
formation are responsible for the observed excess of cosmic rays at ~ 10 18 eV. From 
normalization of the calculated neutron flux to the one observed in the cosmic ray 
excess, we predict the neutrino and gamma-ray fluxes. It has been found that the 1 
km 2 neutrino detector of the IceCube type should detect from a few up to several 
events per year from the Galactic Center, depending on the parameters of the consid¬ 
ered model. Also future systems of Cherenkov telescopes (CANGAROO III, HESS, 
VERITAS) should be able to observe 1-10 TeV 7 -rays from the Galactic Center if 
the pulsar was created inside a huge molecular cloud about 3 — 10 x 10 3 yrs ago. 

Key words: Galaxy: center - Cosmic rays - pulsars: general - ISM: clouds - gamma- 
rays: theory - radiation mechanisms: non-thermal 


1 INTRODUCTION 

Recently the AGASA collaboration has reported the exis¬ 
tence of an extended excess of cosmic rays (CRs) over a 
narrow energy range 10 17 ' 9 — 10 18 ’ 3 eV from directions close 
to the Galactic Center (GC) and the Cygnus region with 
the significance of 4.5<r and 3.9<r, respectively (Hayashida 
et al. 1999). The GC excess was confirmed in the analy¬ 
sis of the SUGAR data (Bellido et al. 2001), in which case 
the observed signal in the energy range 10 1 ' 9 — 10 18,5 eV 
is consistent with that from a point-like source of neutral 
particles, with an estimated flux of (9 ± 3) x 10~ 14 m -2 
s _1 , offset from the true location of the GC by about 7.5°. 
Hayashida et al. (1999) suggested that such an anisotropy of 
CRs in a narrow energy range can be explained naturally by 
neutrons due to the fact that particles with such energies are 
able to reach the Earth before decaying from distance of the 
GC. Numerical simulations of charged particles (protons) 
propagation in reasonable models of the Galactic magnetic 
field give rise to the more extended source which may even 
be significantly shifted in the sky from the original direc¬ 
tion towards the source (e.g. Clay 2000, Bednarek, Giller 
& Zielinska 2001). However, the results of charged hadrons 


propagation strongly depend on the details of the magnetic 
field model which is not well known at present. 

The production of neutrons in discrete galactic sources 
was already discussed some time ago by Jones (1990), Som¬ 
mers & Elbert (1990) and Bednarek (1992) in relation to 
reports of the excess of EeV particles from the direction 
of Cygnus X-3 (Cassiday et al. 1989, Teshima et al. 1990). 
It has recently been suggested that in the GC region neu¬ 
trons might be produced in collisions of hadrons with matter 
(Takahashi & Nagataki 2001). Hadrons can be accelerated 
by a massive black hole associated with the Sgr A* (Levin¬ 
son & Boldt 2000), or by the shock waves of supernovae 
which explode into their own stellar winds (Rhode, Enslin 
& Biermann 1998). 

The Galactic Center region (inner ~ 50 pc) is rich 
in many massive stellar clusters with a few to more than 
100 OB stars (Morris & Serabyn 1996, Blum et al. 2001). 
These stars should soon explode as supernovae. In fact re¬ 
cent multiple supernova explosions in the GC region (10 3 
supernovae in the past 10 s years) are suggested by the ob¬ 
servations of the diffuse hot plasma emitting X-rays (Ya- 
mauchi et al. 1990). For example one remnant of such a 
young supernova with the age of ~ 80 years (GO.570-0.018) 
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has recently been reported by Senda et al. (2001). Since it is 
expected that pulsars are formed in explosions of such mas¬ 
sive stars, we can expect that the GC region should contain 
some young pulsars, a number of them being 10 2 — 10 3 yrs 
old. Motivated by these observational results we may as¬ 
sume that at least one of these young pulsars, formed in a 
supernova explosion of the type Ib/c, has parameters which 
allow the acceleration of iron nuclei to energies ~ 10 2 ° eV, 
as postulated by some models for the Galactic origin of the 
highest energy cosmic rays (Blasi, Epstein & Olinto 2000; 
De Goubeia Dal Pino & Lazarian 2000). We suggest that 
neutrons from desintegration of iron nuclei, which are ac¬ 
celerated by such energetic pulsar in the GC region, can be 
responsible for the observed excess of the cosmic rays with 
energies ~ 10 18 eV. Note that the iron nuclei with energies 
~ 10 2 ° eV and neutrons with energies ~ 2 x 10 18 eV have 
the same Lorentz factors. In order to test this hypothesis 
we predict the neutrino and y-ray fluxes accompaning the 
process of neutron injection by iron nuclei. 


2 A PULSAR INSIDE A MOLECULAR CLOUD 

We investigate the scenario in which a very young pulsar 
is formed in a core collapse of the type Ib/c supernova 
immersed within a huge molecular cloud (or high density 
medium), characteristic of the GC region. Such supernovae 
are probably projenitors of neutron stars with extreme pa¬ 
rameters (period, magnetic held surface), which allow them 
to accelerate iron nuclei to the highest energies observed 
in the CRs. In about a year after supernova explosion, the 
pulsar is surrounded by such a dense expanding envelope 
that hadrons, accelerated in the pulsar magnetosphere or 
the pulsar wind zone, are not injected into the surrounding. 
These hadrons lose energy on multiple inelastic collisions 
producing high energy neutrinos (Protheroe, Bednarek & 
Luo 1998, Bednarek 2001, Beall & Bednarek 2001). In later 
stages the particles can escape into the surrounding and dif¬ 
fuse in the magnetic held of the cloud, suffering collisions 
with the matter from time to time. In our further considera¬ 
tions we discuss, as an example, two media typical of the GC 
region in which the pulsar may be immersed. The hrst one 
is a huge molecular cloud with the radius R c = 10 pc, the 
density n c = 10 3 cm -3 , and the magnetic held B c = 10~ 4 
G ( the total mass ~ 10 5 Mq), and the second one is an 
extended high density region inside the GC with R c = 50 
pc, n c = 10 2 cnh 3 , and B c = 3 x 10 -5 G (the total mass 
~ 10 6 Mq). The clouds with such parameters are typical of 
the GC region and of the Nuclear Bulge, which is a narrow 
layer of the interstellar matter in the central ~ 600 pc of the 
Galaxy with the height of ~ 50 pc. 


2.1 Acceleration of nuclei by a pulsar 

Following the recent works of Blasi, Epstein & Olinto (2000), 
Beall & Bednarek (2001), and Bednarek (2001) we assume 
that pulsars can accelerate iron nuclei in its wind zone 
(the mechanism called magnetic slingshot, Gunn & Os- 
triker 1969). The maximum energies which the nuclei can 
reach in this model, are determined by the magnetic held 
energy per particle in the pulsar wind zone and depend on 


the pulsar parameters as, 

E Fe = . B (r , L0) , « 1.8 x 10 n Bi 2 P“ s 2 GeV, (1) 

87TO G j(r L c) 

where P = 10 _3 P ms s is the pulsar period, B = 10 12 Bi2 G is 
the pulsar’s surface magnetic held, n G j = R(»~Lc)/(2ecP) « 
3.3 x 10 1 4 1?i2 Pins cm 3 is Goldreich & Julian (1969) density, 
r LG = cP/2n « 4.77 x 10 6 P ms cm, and c is the velocity 
of light. According to the slingshot mechanism, the accel¬ 
eration of nuclei occurs very fast so that they do not lose 
energy during this stage. During the hrst ~ 1 yr after explo¬ 
sion these nuclei interact with the radiation held and matter 
of the supernova envelope (Beall & Bednarek 2001). Only 
when the envelope becomes transparent, can the iron nuclei 
be injected into the surrounding of the parent pulsar’s super¬ 
nova. We obtain the spectrum of iron nuclei injected by the 
pulsar following the derivation by Beall & Bednarek (2001). 
Their calculation is based on the general prescription given 
by Blasi, Epstein & Olinto (2000) in which the number of 
nuclei accelerated to energies E scales as a part £ of the Gol¬ 
dreich & Julian (1969) density at the light cylinder radius. 
Beall & Bednarek (2001) modify this approach by noting 
that the spectrum injected by the pulsar at a hxed time 
’t’ should be non-monoenergetic due to the fact that the 
magnetic energy density, responsible for the acceleration of 
particles, changes along the pulsar’s light cylinder height. 
As a result, a pulsar with specific parameters injects nuclei 
with the spectrum which can be below Ep e described by (see 
details in Beall & Bednarek 2001), 


dN _ 27re;r- 3 c n G j(E Fe g 2 )~ 1/3 
iEdt 3[(E Fe /.E) 2 / 3 -l] 1/2 

^ 3xlQ 30 g(B 12 p- a 2 E-() 2 / 3 Fe 

[(B F e/B) 2 / 3 -l] 1/2 s Ge 


( 2 ) 


We estimated that these nuclei can escape through 
the supernova envelope after ~ 1 yr after the supernova 
explosion for typical parameters of the supernova, i.e the 
mass of the envelope in the case of type Ib/c supernovae 
M env = 3 Mq, and the expansion velocity of the envelope 
at the inner radius is v e nv = 3 x 10 s cm s _1 (Beall & 
Bednarek 2001). The iron nuclei diffuse in the magnetic 
held of the high density medium in the GC region, i.e. huge 
molecular clouds. Some of them interact producing neutrons, 
neutrinos, and 7-rays. 

In order to obtain the equilibrium spectrum of iron nu¬ 
clei inside the cloud, we have to integrate over the activ¬ 
ity time of the pulsar since its parameters evolve in time 
due to the pulsar’s energy losses. If we assume that the 
pulsar loses energy only on electromagnetic waves, then 
its period changes according to the formula P^ a (t 0 b s ) = 
1.04 x 10~ 9 t o bs-B?2 + Po,ms, where P 0 ,ms Pms (fobs) are the 
initial and present periods of the pulsar. The equilibrium 
spectrum of iron nuclei at a specific observation time, fobs, 
is calculated from 


dN 

dE 



dN 

dEdt 


Ke- c(t ° bs ~ t)/x dt, 


(3) 


where to — 1 yr, K gives the part of nuclei produced at the 
time’t’ which do not escape from the cloud due to the diffu¬ 
sion and are still present inside the cloud at the time fobs- A 
is the mean free path for collision of the iron nuclei with the 
matter of the cloud. Note that for the parameters consid- 
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Figure 1. The differential spectra of neutrons (dashed curves) 
and 7 -rays (full curves) produced in the interactions of iron nuclei 
with the matter of a molecular cloud with the radius R c = 10 pc, 
density n c = 10 3 cm -3 , and the magnetic field B c = 10“ 4 G, 
at the time after the pulsar formation: t = 10, 10 2 , 10 3 , 3 X 10 3 , 
and 10 4 yr (from the thinnest to the thickest curve). The initial 
parameters of the pulsar are P ms = 2 and B \2 = 6 . 

ered in this paper the value c(t 0 b s — i)/A, in Eq. j|, is always 
less than one. Therefore, a significant part of iron nuclei 
with energies ~ 10 2 ° eV can escape without interaction as 
postulated by the model of Blasi, Epstein & Olinto (2000). 
The value of K is estimated from K = (R c /Ddn ) 3 , where 
D dif = (rLCt/3 ) 1 ^ 2 is the diffusion distance of iron nuclei in 
the magnetic field of the cloud, and 7 ~l is the Larmor radius 
of the iron nuclei with energy E. For the case, Ddif < R c , 
we take K m 1. We assume that the giant molecular cloud 
does not expand. So then the nuclei do not suffer adiabatic 
energy losses. 

2.2 Interaction of nuclei inside the cloud 

The part of iron nuclei confined within the molecular cloud, 
interact with a relatively dense medium suffering desintegra- 
tions and pion energy losses. The pions decay into neutrinos 
and 7 -rays. Applying the equilibrium spectrum of iron nuclei 
(see Eq. pi), we calculate the differential spectra of neutrons 
(from desintegrations of the iron nuclei), muon neutrinos, 
and 7 -rays (from inelastic collisions of iron). In these cal¬ 
culations we assume that: ( 1 ) in a single interaction, the 
iron nuclei desintegrate into two nuclei and the rest of neu¬ 
trons is released with the Lorentz factors comparable to the 
Lorentz factor of the primary nuclei; (2) pions, produced 
in interactions of iron nuclei with the matter of the cloud, 
have Lorentz factors comparable to their parent nuclei. They 
are produced with multiplicities given in Orth & Buffing¬ 
ton (1976). 

The computations of differential neutron, neutrino, and 
7 -ray spectra have been done for the pulsar with the sur¬ 
face magnetic field B 12 = 6 (typical of the observed radio 
pulsars) and the initial period P ms = 2. Such a pulsar is 
able to accelerate nuclei to such high energies (see Eq. |I]) 
that neutrons from their desintegration fulfil the observa¬ 
tional constraint put by the excess of cosmic ray particles 
at 10' 9 — 10 8,5 GeV (Hayashida et al. 1999, Bellido et al. 
2001). The pulsar born with such parameters slows down 
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Figure 2. The differential spectra of neutrons (dashed curves) 
and 7 -rays (full curves) observed from the Galactic Center on 
Earth at the time 10,10 2 ,10 3 ,3 X 10 3 , and 10 4 yrs (from the 
thinnest to the thickest curve) after the formation of the pulsar. 
The initial parameters of the pulsar are as in Fig. Q. The thick 
dotted curve shows schematically the observed cosmic ray spec¬ 
trum (see e.g. Boratav & Watson 2000) within the 20° circle. 

due to the dipole energy losses to 11 ms after 10 2 yr, 34 ms 
(10 3 yr), 60 ms (3 x 10 3 yr), and to 106 ms (10 4 yr). So then 
it is still fast enough to produce pulsaed 7 -rays (Crab and 
Vela type pulsars). 

As an example in Fig. [I] we show the spectra of neutrons 
(dashed curves) and 7 -rays (full curves) produced by the 
iron nuclei inside a cloud with the radius R c = 10 pc, the 
density n c = 10 3 cm -3 (the cloud mass 10 5 Mq), and the 
magnetic field B c = 10 -4 G at different times after the 
pulsar formation tabs = 10 , 10 2 , 10 3 ,3 x 10 3 , and 10 4 yrs. 
As expected, the spectra of neutrons extend up to a few 
10 18 eV. However, their intensities at these highest energies 
significantly drop with time after the pulsar formation due to 
the escape of iron nuclei from the cloud for the time t a b B > 10 
yr, which is caused by their diffusion in the magnetic field. 
Due to the large multiplicities of pion production at these 
high energies, the maximum energies of neutrinos and 7 -rays 
are about 2 orders of magnitude lower than for neutrons. 
Their spectra at the liigiest energies also drop with time 
due to the escape of nuclei from the cloud. 


© RAS, MNRAS 000, 000-000 













4 W. Bednarek 



log (E / GeV) 

Figure 3. The integral spectra of 7 -rays produced in hadronic 
interactions of iron nuclei with the matter of the molecular cloud 
at the Galactic Center for the models considered in the text: I 
(full curve), II (dotted), III (dot-dashed), and IV (dot-dot-dot- 
dashed), and V (dashed). The thick-dashed and full curves show 
the sensitivity of the operating HEGRA System and planned 
CANGAROO III, HESS, and VERITAS systems for a 50-hour 
exposure to a single source (Lorenz 2001, and references therein). 


3 NEUTRONS, NEUTRINOS, AND 
GAMMA-RAYS ON EARTH 

The large distance to the GC, Dgc ~ 8.5 kpc, has influence 
on the expected fluxes of neutrons and 7 -rays measured on 
Earth, due to the decay of unstable neutrons and the ab¬ 
sorption of 7 -rays in the microwave background radiation 
(MBR). The mean free path for neutrons depends only on 
their energy, An = ctn7n ~ 9.8EN,EeV kpc, where a neu¬ 
tron’s lifetime is tn = 918 s, its energy En = 9.38 x 1CU 10 7n 
in EeV, and 7 n is the Lorentz factor. The mean free path for 
the absorption of photons in MBR, A 77 , was computed in 
several papers just after the discovery of MBR (e.g. Gould & 
Schreder 1966). A 77 becomes comparable to the distance of 
the GC at energies of ~ 1 PeV. Therefore the fluxes of neu¬ 
trons and 7 -rays observed on Earth are reduced by the factor 
exp(— Dgc/A 77 ; n), where A 77; n = A 77 or An respectively. In 
Fig. 1 we show the fluxes of neutrons and 7 -rays observed 
on Earth at different times after the pulsar formation for 
the case of the pulsar with the previously mentioned initial 
parameters and two different sets of parameters describing 
the medium in the GC region mentioned in section 2: the 
case of the molecular cloud (the upper figure), and the case 
of a more extended and lower density region around the GC 
(the lower figure). The comparison of Figs. |I| and ^ allows 
us to conclude that the absorption effects of 7 -rays become 
important in the PeV energy range, and that neutrons with 
energies below ~ 10 s GeV are not able to reach the Earth. 
Since the propagation of neutrinos is not influenced by any 
process, their fluxes can be simply obtained for the known 
distance to the GC. In Fig. ^ we mark also the observed 
spectrum of cosmic rays (CR) below 3 x 10 18 eV within the 
20° circle (the analysis box of the AGASA data). We find 
that only the pulsars born within the last ~ 3 x 10 3 yrs 
are able to provide fluxes of neutrons which exceed the CR 
limit, provided that they accelerate nuclei with the efficiency 
£ = 1 (see Eq. |). 


Table 1. Gamma-rays and neutrinos from the Galactic Center. 


Model 

N 7 (> 1 TeV) 

iV 7 (> 10 TeV) 

N a 

1 y V 

AT na 

1 y V 

(I) 

4.3 x 10 “ 12 

2.2 x 10~ 12 

23 

30 

(II) 

8.7 x 10 “ 13 

6.6 x 10- 13 

11 

16 

(III) 

1.4 x 10 “ 13 

1.25 x 10 “ 13 

4.2 

7.1 

(IV) 

2.5 x 10 -13 

1.7 x 10 -13 

5.3 

0 ° 

00 

(V) 

6.7 x 10 “ 14 

5.5 x 10 ~ 14 

2.0 

3.8 


4 DISCUSSION AND CONCLUSION 

The SUGAR group estimates the flux of particles which 
causes reported excess of the cosmic ray particles in the 
energy range 10 17 ' 9 — 10 18,5 eV on (9 ± 3) x 10 ~ 14 m -3 s _1 
(Bellido et al. 2001). If this excess is caused by neutrons pro¬ 
duced in the pulsar model discussed here, then the expected 
flux of neutrons can be compared with the observed one. 
Basing on this normalization we predict the fluxes of neu¬ 
trinos and gamma-rays on Earth. This procedure allows us 
to derive the free parameter of our model (i.e. the efficiency 
£ of iron acceleration by the pulsar) and limit the age of the 
pulsar for other fixed parameters, P, B, R c , n c , B c , which are 
in fact constrained by the observations. We consider five dif¬ 
ferent sets of parameters describing our scenario: model (I) 
R = 10 pc, n = 10 3 cm -3 , B c = 10 -4 G, t 0 b s = 10 4 yr; 

(II) fobs = 3 x 10 3 yr and other parameters as above; (III) 
lobs = 10 3 yr and other parameters as above; (IV) R = 50 
pc, n = 10 2 cm -3 , B c = 3 x 1CU 5 G, t 0 b s = 10 4 yr; and 
(V) Cbs = 3 x 10 3 yr and other parameters as in (IV). They 
all concern two sets of parameters for the medium in which 
the pulsar is formed, and differ in the pulsar’s age which is 
not constrained by any observations. In all these models we 
assume that the pulsar is born with B = 6 x 10 12 G and 
P 0 = 2 ms. Normalizing the predicted neutron flux to the 
observed excess of CR particles we derive the value of the 
parameter £ which has to be £ ~ 1 (model I), 0.18 (II), 0.03 

(III) , 0.3 (IV), and 0.09 (V). 

Using the above estimates for £ we can now predict the 
expected fluxes of 7 -rays and muon neutrinos and antineu¬ 
trinos in the case of every model. The integral spectra of 
7 -rays from the GC region are presented in Fig. together 
with the sensitivities of the present HEGRA telescope sys¬ 
tem and the planned next generation of telescopes, i.e. CAN¬ 
GAROO III, HESS, VERITAS. We also report in Table § 
the 7 -ray fluxes above 1 TeV and 10 TeV in units cm~ 
s _1 . Although the 7 -ray spectra have a maximum above 10 
TeV for all models, the 7 -ray fluxes in the energy range 1- 
10 TeV produced in models, (I) ~ 2.1 x 1CP 12 cm -2 s _1 , 
and (II) ~ 2.1 x 10 -13 cm -2 s _ , and probably also in 

(IV) ~ 8 x 10 -14 cm -2 s~\ should be observed by the fu¬ 
ture systems of Cherenkov telescopes of the CANGAROO 
III, HESS, and VERITAS type. Models (III) and (V) pre¬ 
dict fluxes below the sensitivity limit of these Observatories. 
Only the HEGRA Collaboration observed the Galactic disk 
including the GC region (Pohlhofer et al. 1999). The upper 
limit on the possible sources, equals 1/4 Crab in the Galac¬ 
tic plane which is above the 7 -ray flux predicted even by the 
model (I). However, since the GC region can be observed by 
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Figure 4. The differential spectra of muon neutrinos and an¬ 
tineutrinos produced in hadronic interactions of iron nuclei with 
the matter of a molecular cloud at the Galactic Center for the 
models discussed in Fig. 3. The dashed curves indicate the atmo¬ 
spheric neutrino background, ANB, (Lipari 1993) within a 1° of 
the source and the dotted line shows the 3 yr sensitivity of the 
IceCube detector (Hill 2001). 


this experiment only at zenith angles larger than 60°, this 
limit does not refer to the GC region. 

In Fig. |I]we show the muon neutrino and anti-neutrino 
spectra expected in the discussed model. At energies > 10 
TeV, these spectra are above the expected flux of atmo¬ 
spheric neutrino background (ANB) and also above the 3 yr 
sensitivity limit of the planned large size neutrino detector 
IceCube (Hill 2001). We estimate the number of muon neu¬ 
trino events during one year in the IceCube detector basing 
on the calculations of the likelihood of detecting such neu¬ 
trinos by a detector with a surface area of 1 km 2 obtained 
by Gaisser & Grillo (1987). The results of these calculations 
are shown in Table 1. We distinguish the case of neutrinos 
coming to the neutrino detector from directions close to the 
horizon, i.e. not absorbed by the Earth ( N “), and neutri¬ 
nos which arrive moving upward from the nadir direction, 
i.e. partially absorbed by the Earth (A®) (for absorption co¬ 
efficients see Gandhi 2000). From Table 1 it is clear that the 
IceCube detector should detect a few up to several neutri¬ 
nos per year from the Galactic Center region provided that 
the excess of cosmic rays at ~ 10 18 eV from the GC region 
is caused by neutrons from desintegrations of iron nuclei, 
accelerated by a very fast pulsar. The detection of the pre¬ 
dicted fluxes of neutrinos from the Galactic Center (or lack) 
will also place constraints on the recent model of extremely 
high energy cosmic ray production in the pulsar scenario 
(Blasi, Epstein & Olinto 2000), since the parent iron nuclei 
which inject neutrons with energies ~ 10 18 eV, have to be 
accelerated to energies ~ 10 20 eV. 

The EGRET detector on board the Compton GRO has 
detected a strong 7 -ray source with luminosity of ~ 2 x 10 3 ' 
erg s _1 (Mayer-Hasselwander et al. 1998). This emission 
seems to come from the extended region with the radius 
of ~ 80 pc around the GC. In another paper (in prepara¬ 
tion) we suggest that this emission can be explained in terms 
of the general scenario discussed here if a part of energy of 
the accelerated iron nuclei is transferred to the relativis¬ 
tic positrons due to the resonant scattering in the pulsar 


shock region (see Gallant & Arons 1994). These positrons, 
accumulated in the cloud, produce high energy radiation in 
synchrotron and inverse Compton processes. 
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